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(54) A method and system for detecting grids In a digital image 

(57) A method for detecting grids in R digital image 
havings plurality of pixels tor capturing an image, com- 
prises" the steps of performing a high pass tiller opera- 
lion on a sample of pixels lor filtering out the image from 
the sampled pixels; and performing a statistical F-tesi 
operation on the results of the high pass filter operation 
for determining when the digital image conlains grids. 



□ 



f "" i 

s 


1 — 

p 













\ \ \ i - / 

\\T1 ; 



. 1 r 



2\ 



4ir 



rnoiedcj jouv* 7z>x>\ --Anio <rn, 



DEC 13 2001 13:34 



PAGE. 53 



EP 0 723 782 A2 



2 



Description 

Background Of The invention 

Field of invention 

The invention relates generally io a method and 
system for detecting grids in a digital image. 

Background of the Related Art 

There are methods and systems, which are well 
known in the art. for viewing internal portions cf an object 
by directing radiation signals, generally referred to in the 
art as x-rays, onto the object and by receiving the radi- 
ation signals which pass through the object onto a re- 
cording medium, typically film, ideally, the radiation sig- 
nals will pass from the radiation source to the mm in a 
straight line. However, as the radiation signals pass 
through the object, some of the s»gnyis are scattered so 
that they no longer tiavel in a stiaight line, which scat- 
tering causes reduced image quality. To reduce this ef- 
fect, radiographic grids are placed between the object 
cjnd film. Radiation grids are lead foil strips placed in a 
spaced apart relationship with each other, which spaced 
apart relationship is typically maintained by placing alu- 
minum between the strips. The strips and aluminum are 
then packaged in a suitable covering. These grids, 
though effective for reducing rad.ation scattering, occa- 
sionally introduce artifacts such as grid lines into the im- 
age. 

The film containing the captured image is proc- 
essed in a variety of well known and utilized techniques 
for permitting on observer to view the image captured 
on the film. One such technique is to convert the image 
into a digitized form for displaying the image on a com- 
puter display and similar devices 

To convert the image into digitized form, typically, a 
digitizer converts the film into digital form for permitting 
the image to be displayed on a computer display by well 
known techniques. 

Before displaying the image on the computer 
screen, grids lines, if there are any. shculdbe removed. 
There are two well known mcthocs tn the art tor reducing 
tne eftects of grid lines. One methoc utilizes moving the 
grids when the object of interest is radiated. The other 
method is to apply a processing algorithm to the com- 
pute! system which is displaying the image. 

Although the above methods foi reducing gud lines 
are satisfactory, they are not without drawbacks. Moving 
grids involve a costly procedure and, in addition, require 
the object to be exposed to the radiation for a longer 
period of lime than iG normally required. This obviously 
is a drawback if the object is a person. For these rea- 
sons they arc not a preferred method of reducing the 
effects of grid lines. The processing algorithm is also 
costly and. as a secondary effect can reduce the quality 
of the image 



Consequently a need exists for a method o» detect- 
ing grids so that costly processing algorithms for remov- 
ing grids are applied only when there are grids in the 
image. 

s 

Summary Of The Invention 

It is an object of the present invention to provide a 
method and system for detecting grids 

70 u is an advantage of th e present invention to provide 
such a grid detection method and system which are cost 
effective and involve a minimum of time to implement. 

It is a feature of the p*esent invention to implement 
a high pass filter for removing the capture image from 

is the digitized representation and to apply a statistical op- 
eration on me results for detecting the presence of grids. 

The above and other objects, features and advan- 
tages ot the present invention will become more appar- 
ent when taken in conjunction with tne Allowing descrip- 

20 Uon and drawings wherein like characters indicate like 
paits and which drawings form a part of the present in- 
vention. 

With this object in view, the present invention re- 
sides in a method for detecting grids m a digital image 

55 having a plurality of pixels for capturing a^ image, com- 
prising the steps of (a) performing a high pass filter op- 
eration on a sample of pixels for filtering out the image 
from the sampled pixels: 2nd (b) performing a predeter- 
mined mathematical operation on the results of the high 

30 pass filter operation for determining when the digital im- 
age contains grids. 

Brief Description Of The Drawings 

55 in the course of the following detailed description, 
reference will be made to the anached drawings in 
which: 

Fig. 1 is a drawing illustrating the environment of 
the present invention 

Fig. 2 is a flowchart of the method in accordance 
with the present invention. 

Fig. 3 is an example of a plot pixel value versus pixel 
location: 

^ Fig. 4 is a plot in which Fig. 3 is smoothed: and 

Fig. 5 is tne subtraction of Fig. 4 from Fig 3. a high 
pass fine: 

Detailed Descripion Of The Invention 

so 

Referring to Fig. 1, there is illustrated a system for 
implementing the present invention. The system in- 
cludes a high vo tagc generator 10 which powers an X- 
ray tube 20 for generating a beam of radiation signals 
53 30. generally referred to as X-rays for irradiating an ob- 
ject 40. The radiation signals pass through the objeel 
and are detected by a recording medium, typically film 
50 As is illustrated by tne dashed lines, some Of the 
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radiauon signals are scattered as they pass through the 
obied <0 These scattered signals reduce IN. i image 
quality, and as a result, a well known and uMaod radio- 
graphic qrid 60 is placed between the object 40 and him 
50 lor reducing the effects of scattering. 

The film 50. which now contains the captured im- 
age, is Placed in a digitizer 70 lor convening" into a 
digital formal During this conversion, the digger 70 
2ns one portion el the film 50. generally referred l to m 
,he art as a pxel 80. and assign* a pixel value , to he 
pixel 80. which value represents the contents on the wn 
50 at that particular pixel 80, The digitizer 70 scans the 
entire image in the above described manner tor conven- 
ing each portion of the film 50. or pixel 60 into a pixel 
value In essence, the digitized data represents each 
pixel SO location and its corresponding pixel value. The 
Mm 50 is typically scanned in such a manner that m pixel 
reoresentation the image includes a plurality ot contig- 
uous pixels 80 having a checkered array configuration. 

The dtgtti/ed data is passed lo an im*ge processor 
90 wheie the data is stoied in pixel lotma'. in a f.ist iQO 
or second memory 110. In other words, a pixel location 
and i:s associated pixel value arc stored in either ol the 
two memories 100 and 110. The image processor 90 
also cor.ta.ns an arithmetic unit 120 lor pcrlorming arith- 
m»tie f uncfions and other typical software (unctions. The 
,mage processor 90 upon command from a user can dis- 
play the image stored in the memories 1 00 and 110 on 
a computet display 130 

However, belore the image is displayed on the 
screen 130. a method ol the present invention analyzes 
the image lor the presence of grids which may have 
been created by the radiographic grids 60. Relemng lo 
Fig 2 the method is illustrated in flowchart form. The 
image processor 90 loads a software program ol the 
present invention wh:ch randomly samples 1*0 two per- 
cent ol the oixel values stored in memory 100 and HO 
for determining the mean value 1 50 of the sampled pix- 
els 80. which mean value is an estimate Ol the mean Of 
the entire, image it is instructive to note that two percent 
sampling is mere y the percentage ol the preferred em- 
bodiment and that any percentage of the pixels 80 may. 
in fact be sampled for this estimation. Pixels 80 with a 
value ol me mean or higher are an est.mat.on ol where 
.1 is easier to detect grids Next, an array ol pixels 80. 
preferably a 16 x 16 array, having a pixe value ot the 
mean or higher is selected 160 for determining the mean 
vslue 170 ol the 16 x 16 anay. Again, although a 16 x 
15 anay is preferred, any size anay may be used for 
this sample. If the array mean is below the estimated 
mean ol the image 180, the array is discarded and an- 
other 15 x 16 array is sampled 160. If the array mean .s 
at or above the estimated mean of the entire image loO. 
the array will be further analyzed tor the presence of 
qrids. 

~ The 16 x 16 array is further sampled by sampling 
two or more lines of pixels which are parallel to an im- 
aginary x-axisor ordinate 190 which is through the cent- 



er ol the array. Each pixel 80 from each of the sampled 
„nes whieh h»ve a same ordinate is then averaged 200 
together so that the result is an average pixel value 
atong the ordinate. This result is plotted as a graph ol 
5 pixel value versus p.xel location along the or*n* An 
example o. such a plot is shown m Fig. 3. Referring back 
" Z flowchart. th» Plot is smoothed 210 wh.ch » we" 
KnLvvn n tne art and an illustration ol which ,s shown in 
STT (See DWalPicwre Processing Picture. Bosen- 
,o fe d and Kak. Academic Press ig82 lor further dsclo- 
sure ol smoothing) The smoothed plot is sub.rac.ed 220 
from the ave.aqed plot lor forming a high pass Uter of 
.h^esultant pixels. An example of a plot ot the h.gh pass 
filter is shown in Fig 5. 
,5 Relerrinc to Fiq. 2. the resultant high pass plot is 
autocorre.ate"d 230 by tcchmques which are well known 
fntheart.orenhanc : ng.heh,ghpassplot(SeeD e! ec^rt 

cf Signils in Nc*«. Whaien 1971 for lunher disclosure 
of autocorrelat.on) For further enhancement, the auto- 
,o correlated 240 plot is again aulocorre.a.ed for unhe 
still bene, quatitv. The va,.ance which ,s a well known 
Statistical technique is then calculated 230 from the re- 
-St ot this autocorrelation (See ProbeMty and Staus- 
Deeroet Add-son Wesley 1975 lo, details o. van- 
25 aocc calculation) 

Steps 190 through 250 are then repealed lor sam- 
p,ing along the y-axis or absc ssa. and are recited below 
Tor thoroughness oi understanding. F***^J%?? 
x 1 6 array further samples two or lines of p.xels which 
so are pai to an imaginary the abscissa 260 «h«h» 
also through the center of the array. Each p.xel from 
each of the sampled lines which have a same abscissa 
is then averaged 270 together so that the result • an 
average pixel value along the abscissa This esu t is 
ss plotted as a graph ol pixel value versus pixel locator, 
alono the abscissa and the plot is smoothed 2*0-The 
smoothed plci is subtracted 290 from the averaged plot 
for lormmg a h.gh pass W.er o. the resultant pixels. 
The resultant high pass plot is autocorrelaiec 300 
40 for enhancing the high pass plot. For ^^JJS 
ment. the autccorreiated plot .s aga.n auioccrrelated 
310 for lunher still better quality. The variance is then 
calculated 320 from the result ol this autocorrelation 
After the vanances m both directions are calculated. 
as a F-iest 330 .s used tor determining whether g"ds are 
present in the 16 x 16 sampled array Of pixels ,he F- 
test is calculated as foiiov/s. 



r.'Jjo riiontal 



F= i^jj^.!i^J^ Equation!. 1) 

o 2 Vertical 

SO 

If the F-test score is greater than 1 this results in a de- 
termination that vertical grids are present in xhm »"*P'«° 
array, and a vote ol positive one is given, if the F-tc- 
score is less than 1 this .esuits in a determinate that 
ss horizontal grids are p.esent ,n the sampled array and a 
vote Of negative one is grven. A F-test score c zero 
means that nei'her vertical nor hon/ontal grids are 
present, anc; a vote of 7ero is giver, 
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if all of a pf ©determined percentage of pixels 30 with 
a mean value at or greater than me mean pixel value of 
the entire image have been included in a 16 x IS array 
for analysis 340. the average of all the votes 350 are 
taken and the result indicates whether any grids are 
present and. if so. what direction they rest. The result of 
the average has the same meaning as the F-test score: 
posilve one means vertical grids, negative one means 
horizontal grids are present, and a zero indicates no 
grids are present, if ell the pixels have not been included 
in a 16 x 16 array for analysis 340, the above steps are 
repeated until they have been included in a 1 6 x 1 6 array 
analysis. 

It is also obvious for the man skilled in the art that 
said inventton could also be used lor circular grids, in 
such a case the average of the pixel value is made ao 
cording to a predetermined pattern dependent on the 
location of the 16x16 sampled array with respect to the 
cenier of the grid 

Referring back to Fig. 1, if a grid is present in the 
image, the image processoi 90 applies a filter to the im- 
age for removing them. SJCh filters are well known in 
the art. 



Claims 

1. A method for detecting grids in a X-Ray digital im- 
age having a plurality of pixels for capturing an im- 
age, comprising the steps of 

(a) performing a high pass filter operation on a 
sample of pixels for filtering out the image from 
the sampled pixels: and 

(b) performing a predetermined mathematical 
operation on the results of the high pass filter 
operation for determining when the digital im- 
age contains grids 

2. The method as in claim 1 wherein the predeter- 
mined mathematical operation includes a F-test 

3. The method as in claim 1 wherein the high pass filter 
operation includes the steps of: 

(ai) forming a plot of pixel location versus the 
pixel value averaged, at the corresponding pix- 
el location, according to a predetermined pal- 
tern dependant on the ujid used. 
(32) smoothing the results of the plotting oper- 
ation; and 

(a3) subtracting each plotted pixel value ©nd 
each smoothed pixel value which arc from the 
same pixel location for obtaining the high pass 
filter. 

4. The method as in claim 3 wherein step (al) includes: 



6 

(ai) further sampling from the sampled pixels at 
least two spaced ap^rt lines of pixels: 
(aii) averaging pixel values from each of the 
spaced apart lines which have a matching ab- 
5 scissa for removing noise from the digital im- 

age: and 

(aiii) forming the plot in step (al) from the aver- 
aged pixel values. 

to 5. The method as in claim S wherein step (aii) includes 
averaging pixel values from each of the spaced 
apart lines which have a matching ordinate for re- 
moving noise from the digital image. 

'5 6. The method as in claim 5 wherein the plot is auto- 
correlated for enhancing the plot. 

7. The method as in claim 2 wherein the F-test is F= 
horizontal variance/vertical variance. 

20 

8. The method as in claim 7 wherein vertical grid lines 
result in the F-test having the horizontal variance 
less than the vertical variance, and horizontal grid 
lines have the horizontal variance greater than the 

55 vertical variance. 

9. The method as in claim 3 wherein no grids are 
present in the digital image when either of the con- 
ditions in step & are not met. 



90 

10. The method as in claim 9 wherein each of the sam- 
pled pixels designate a first signal if vertical grids 
are detected and designate a second signal if hor- 
izontal grids are detected for determining which of 
ss me first or second signals have a majority which, in 
turn, determines whether the digi;al image has ver- 
tical or horizontal signals. 
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(54) A method and system for detecting grids in a digital image 



(57) A method lor detecting grids in a digital image 
having e plurality of pixels tor capturing an image, com- 
prises the steps of performing a hish pass filter opera- 



tion on a sample of pixels -or filtering out the imagejrom 
the sampled pixels: and performing a statistical F-test 
operation on :he results o 2 the high pass filter operation 
for determining when the digital image contains grids. 
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